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Abstract 

Enlightened by the idea of the 3x3 CKM angle matrix proposed recently by Harrison et al, 
we introduce the Dirac angle matrix <I> and the Majorana angle matrix ^ in the lepton sector for 
Dirac and Majorana neutrinos respectively. We show that in the presence of CP violation, the angle 
matrix $ or ^' is entirely equivalent to the complex MNS matrix V itself, but has the advantage of 
being real, phase rephasing invariant, directly associated to the leptonic unitarity triangles (UTs) 
and do not depend on any particular parametrization of V. In this paper, we further analyzed how 
the angle matrices evolve with the energy scale. The one-loop Renormalization Group Equations 
(RGEs) of <I>, ^' and some other rephasing invariant parameters are derived and the numerical 
analysis is performed to compare between the case of Dirac and Majorana neutrinos. Different 
neutrino mass spectra are taken into account in our calculation. We find that apparently different 
from the case of Dirac neutrinos, for Majorana neutrinos the RG-evolutions of ^' and strongly 
depend on the Majorana-type CP-violating parameters and are more sensitive to the sign of Am|^. 
They may receive significant radiative corrections in the MSSM with large tan /3 if three neutrino 
masses are nearly degenerate. 
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I. INTRODUCTION 



Since 1998, a number of successful neutrino oscillation experiments have provided us with 
very convincing evidence that neutrinos are massive and lepton flavors are mixed The 
flavor mixing among three neutrinos can be described by the Maki-Nakagawa-Sakata (MNS) 
matrix V. Now three mixing angles in V and two squared neutrino mass differences have 
been approximately determined. But whether neutrinos are Dirac or Majorana particles 
remains an open question. The Majorana nature of massive neutrinos can be revealed 
by the investigation of processes in which the total lepton charge L changes by two units 
AL = 2. Neutrinoless double-beta (Oi//3/3) decay experiments are considered as the most 
promising method in this catalog. If Oz//3/3 decay is eventually observed, we shall make sure 
that neutrinos are Majorana particles. If there is no experimental signal for the Oi//3/3 decay, 
however, we shall be unable to conclude that neutrinos are just Dirac particles 

CP violation in the lepton sector is another open question. In the framework of three 
Dirac neutrinos, CP violation in the MNS matrix V can be described by a single Dirac 
CP-violating phase which can be measured in the neutrino oscillation experiments. If three 
neutrinos are Majorana particles, two extra Majorana CP-violating phases are introduced in 
V. It is well known that the presence of Majorana phases introduces some novel features in 
leptonic CP violation, like the possibility of having CP violation in the case of two Majorana 
neutrinos 3| as well as having CP breaking even in the limit of three exactly degenerate 
^. Besides, these extra Majorana phases can affect significantly the rates of 

ay an important role in the RG-evolutions of 
7|, and be the source of CP violation in the 



neutrmos 



Ou/S/S decay jst and some LFV decays {gI, p^ 



the neutrino masses and mixing matrix V 

n 

leptogenesis [8]. Although the Majorana phases can not be directly measured in neutrino 
oscillation experiments [9], the constrains on them can be drawn indirectly from the studies 
of above mentioned processes. 

If the 3x3 MNS matrix V is unitary, its nine elements satisfy the following normalization 
and orthogonality conditions: 

a i 

where the Greek and Latin subscripts run over {e,fi,T) and (1,2,3), respectively. The six 
orthogonality relations geometrically define six leptonic unitarity triangles (UTs) in the 
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complex plane 10|], as illustrated in Fig. 1. The six UTs have eighteen different sides and 
nine different inner angles, but their areas are all identical to j7/2 with J being the Jarlskog 



invariant of CP violation 



11| defined by 



7 k 

The leptonic UTs provide us with a visible way of studying the CP violation in the lepton 
sector. In the presence of CP violation, all the physical parameters in V can be drawn from 



these six UTs no matter if the neutrinos are Dirac or Majorana particles 



Id, 



18|- 



Enlightened by the idea of the 3x3 CKM angle matrix proposed recently by Harrison et 



al. 



12| . we introduce the 3x3 Dirac angle matrix $ in the lepton sector for Dirac neutrinos. 



Moreover, we extrapolate this concept to account for Majorana neutrinos and propose the 
Majorana angle matrix \E'. In section II, we introduce the concepts of $ and \E'. We show 
that in the presence of CP violation, the angle matrix $ {^) is entirely equivalent to the 
complex mixing matrix V itself for Dirac (Majorana) neutrinos, but has the advantage of 
being real, phase rephasing invariant, directly associated to the leptonic UTs and do not 
depend on any particular parametrization of the MNS matrix. In Section III, we further 
analyze how the angle matrices and some other rephasing invariant parameters evolve with 
the energy scale. The one-loop Renormalization Group Equations (RGEs) of \V\, $, ^ and 
the Jarlskog J are derived. Unlike the CKM angle matrix which is quite stable against 
the RG-evolution [isl, the leptonic angle matrix $ and \E' may receive significant radiative 
corrections when evolve from the electroweak scale to a superhigh energy scale. Section 
IV is devoted to a numerical analysis of the RGE running behaviours of $ and \E', and to 
a careful comparison between the case of Dirac and Majorana neutrinos. Different neutrino 
mass spectra are taken into account in our calculation. A brief summary of the main results 
is given in section V. 



II. DIRAC ANGLE MATRIX V.S. MAJORANA ANGLE MATRIX 



In the mass eigenstate basis, the charged-current interaction of leptons is described by 

9 
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V^3/ 



W- + h.c. 



(3) 



3 



and the CP violation is naturally included if the MNS matrix V is complex. If the neutrinos 
are Dirac particles, one has the freedom to make phase rotations on both the charged lepton 
fields and the neutrino fields, which leads to the redefinition 

V e^'^^Ve-''^-^ , (4) 

where = diag(</);^, </);^, </);^) and = diag((f)j^^, (p^^, cf)^^^. Physical quantities are basis 
independent and must be invariant under the rephasing of Eq. (4). The simplest rephasing 
invariant quantities are the moduli of nine elements of V which are all real. Although 
they might implicitly involve CP-violating phases, it is convenient to look for imaginary 
parameters that explicitly require CP violation. As been pointed by many references, the 
lowest-order (in V) rephasing invariants that are not automatically real are the quartic 
products (which are also called "boxes" in some references) [l^ : 

^ v^,v^^v:^v;, , (5) 

where the Greek and Latin subscripts run over {e,fi,r) and (1,2,3), respectively. are 
not automatically real if a ^ (3 and i ^ j- In this paper, no summation on repeated indices 
is implied. The imaginary parts of "^D/jj 

"^53^,. ^ Im-D^^. ^ Im [V^,V^^V:^Vp\] , (6) 

are the measures of the CP violation in V. For three fermion generations, there is only one 
such independent complex quantity, all ""^^pj equal to the Jarlskog invariant JT" except a sign 
difference. Correspondingly we define the real parts of "'D^^- as 

"^3fJ,,. ^ Re-^D^^. ^ Re [V,,VfSjV:,V^*,] . (7) 

Neutrino oscillation probabilities are linear in "*3?^j and enabling a straightforward 

description of oscillation data. 

j>i j>i 

= 5«,-4E"^5^/..sin^^ 

j>i 

e.,, sin sin sin . (8) 
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General analyses of neutrino oscillations among three flavors can readily determine the boxes 



15| 
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^e2 






Ml 








i^-^rl 


^r2 





Enlightened by the idea of the 3x3 CKM angle matrix proposed by Harrison et al. 12| . 
we can similarly construct the 3x3 MNS angle matrix $ for Dirac neutrinos 



(9) 



by using the angles of nine box invariants 

^ -arg(-^^n^,) , (10) 

where a, /5 and 7 run co-cyclically over e, fi and r, while i, j and k run co-cyclically over 1, 
2 and 3. 

We can easily find that the absolute value of is just the inner angle shared by the 
UTs Aq, and as one can see in Fig. 1. Besides, the common sign of all nine angle matrix 
elements is just the sign of the Jarlskog invariant J^. \i J > 0, all the nine angles in $ lie 
between and vr; if ^7 < 0, all the nine angles lie between — vr and ^. Each row or column 
of $ corresponds to one UT, and the unitarity of V now implies that elements of $ satisfy 
the normalization conditions 

E'^- = E'^- = ^- (11) 

a i 

We can draw from Eq. (11) that there are only four independent real parameters in $, same 
as the number of the independent real parameters in the unitary MNS matrix V . We can 
further prove that in the presence of CP violation, the angle matrix is fully equivalent to 
the MNS matrix. In Appendix A, we show how to re-obtain the mixing matrix V from the 



angle matrix $, the process is analogy to that in the quark sector 



12| 



^ There is another immediate way to find out the sign of J . We firstly define the sequence of the sides 
of the UTs: for triangles Ag, and A^, we follow the sequence of V^iV^^ ^02^/32 ~^ Ki3^^*3) 
where a/3 = /ir, re and e/i respectively; and for triangles Aj^, Aj and A3, we follow the sequence of 
^ei^ej ~^ ^tii^ilj ~^ ^Tt^Tjj where ij — 23, 31 and 12 respectively. Then all the UTs can be sorted into 
two classes: clockwise triangles and anti-clockwise triangles. We can easily find that if > 0, all six UTs 
(as shown in Fig. 1) are clockwise triangles while a negative indicates that all six UTs are anti-clockwise 
triangles. This rule is true for both the Dirac and the Majorana neutrinos. 
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The latest global analysis of current neutrino oscillation data yields 0.27 < sin^ < 0.36, 
0.39 < sin2^23 < 0.64 and 0.001 < sin^^^g < 0.035 (NH) or 0.001 < sin^^^g < 0.039 (IH) 
at the 3a level [l^, where "NH" and "IH" correspond to the normal and inverted neutrino 
mass hierarchies respectively. The CP-violating phases remain totally unconstrained. Cor- 
respondingly, the allowed range of the moduli of the elements of the MNS matrix and 
the elements of the Dirac angle matrix ^^.^ can be obtained: 
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(13) 



where the upper (lower) row corresponds to normal (inverted) neutrino mass hierarchy. The 
Jarlskog invariant J can range between [-0.0433, 0.0433] (NH) or [-0.0455, 0.0455] (IH) at 
the 3a level. 

The question is more pressing when we consider the case of Majorana neutrinos. Due to 
the Majorana nature of the neutrinos, the phases of three neutrino fields in Eq. (3) can not 
be freely chosen. The phase rotations on the charged lepton fields lead to the redefinition 



(14) 



where = diag(0;^, 
invariants llj 



). Therefore, for Majorana neutrinos, we have the new rephasing 

(15) 



c — y Y* 



which are not quartic but quadric products of V^j and are not obviously real if i 7^ j. 
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In order to include the informations of Majorana-type CP violation, we introduce the 
following 3x3 Major ana angle matrix 



(16) 





el 


^e2 






Ml 










^.2 





with its nine elements are the angles of S parameters in Eq. (15) 



(17) 



where i, j and k run co-cyclically over 1, 2 and 3. We can find that the three matrix elements 
in each row of \l/ sum to zero: 



/ ^ cti 



arg 



. 



(18) 



Above normalization conditions are satisfied independent of if V is unitary. Then the number 
of independent real parameters in \E' is six, equals to the number of free parameters in the 
Majorana neutrino mixing matrix. In the case of Majorana neutrinos, we can also reconstruct 
the leptonic mixing matrix V from the Majorana angle matrix \E', the detail processes can 
be found in Appendix A. 

We can easily see from Eqs. (5) and (15) that, there are the constitutive relations 



ai\ — I Q Q* 

— '^aij'^fiij 



(19) 



Therefore we can easily construct the $-matrix from the \l/-matrix: 

$ ■ = \l/ ■ — ^! o- + Tl 



(20) 



where a, /3 and 7 run co-cyclically over e, and r. It means we can also write out the 
Dirac angle matrix $ for Majorana neutrinos, but are unable to draw informations about 
the Majorana phases from it. 

From Eq. (19) we can find that even all "^Q'^j are zero (i.e., JT" = and no Dirac-type 
CP violation in V), the imaginary part of some S^^^ can be nonzero and stands for the CP 
violation in V . On the contrary, if S^^^ are all real, there is no CP violation in V and J 
also equals to zero. It means that Dirac-type CP violation requires the existence of the 



Majorana-type CP violation but obviously the converse is not true 
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Here we introduce another combination of S„,, 

Val3ij = ^aij^Pij = '^aiVaj^pypj ' (21) 

in which informations of the Majorana-type CP violation are also involved. For the sake of 
concision, we define the following notations 

^apij = Re ri^f^ij and I^^^^- = Im r/„^.^. . (22) 

The parameters S^^j and ?7Q,/3ij "^^^ show up in a variety of lepton number violating 
processes including Oi//3/3 decay and possibly leptogenesis et al. For example, the effective 
neutrino mass {m)^^ in Oz//3/3 decay can be expressed as 



iKilV^'^-m, + iK/e-^^*-^^ + iKsr^sl • (23) 



Another point worth to mention is that the nine elements of \& are physically related 
to the orientations of the nine sides of UTs A^^, and A3 in the complex plane (see 
Fig. 1). Therefore, the orientations of these three UTs have physical meanings if neutrinos 



are Majorana particles [18|. On the other side, if the neutrinos are Dirac particles, the 
orientations of the UTs have no physical meaning, reflecting the fact that Dirac UTs rotate 
under rephasing of the charged lepton fields or the neutrino fields. 

Before our analysis of the radiative corrections to the angle matrix $ and we show 
here that the conceptions of the Dirac and the Majorana angle matrices can be extrapolated 
to account for any generation of neutrinos. There is no a unique way to choose the elements 
of $(^) and for (A^ > 3) generation neutrinos, we give one possible choice here. Here 
(e, /X, r, S4, S5, . . . , s^v) stand for the flavor indices and (1, 2, 3, 4, 5, ... , A^) for the mass 
indices. 
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For N generation of Dirac neutrinos, we define the N x N Dirac angle matrix 
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_arg(-'^4Q^_,) 

- arg (-^^□,,5) 
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V - arg (-^^D^a) - arg (-^=^n^4) - arg (-^"D^s) 



and its elements satisfy the normalization conditions 



— arg 

-arg(--in,^2) 
-arg(--^in,^2) 

-argpin^2)y 



(24) 



(25) 



where a stands for the flavor index and i for the mass index. Then there are altogether 
- (A^ — 1) (A^ — 2) independent parameters in which is equivalent to the number of the 
independent real parameters in a A'^ x A'^ unitary mixing matrix. 

For N generation of Majorana neutrinos, the N x N Majorana angle matrix can be 
defined as 
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arg (^4^*5) 
arg {V^,V;,) 
arg {V,,V;,) 



arg(V;iK*2) \ 
arg {V^,V;,) 

arg (K1K2) 



\arg(K.2K;3) arg(K^3F;J arg (F^^.^^s) ••• ^^s{Vs^iV:JJ 



(26) 
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and the matrix elements in each row of satisfy the normahzation conditions 




where a stands for the flavor index and i for the mass index. There are totally -A^ (A^ — 1) 
independent real parameters in Again, this number is equivalent to that of the N x N 

Majorana neutrino mixing matrix. In case of A^ generation Majorana neutrinos, we can also 
construct the $ matrix from the \E' matrix: 

= + (28) 

where /3 and 7 are the next two flavor indices right after a. 



III. ONE-LOOP RENORMALIZATION GROUP EQUATIONS 



The RGEs of neutrino masses and mixing have been discussed in many papers with 
variety of parametrizations [3, [lol-Sl, 23]. It has been shown that three neutrino masses 
and the mixing matrix may receive large radiative corrections, especially if neutrino masses 
are nearly degenerate or in case of the MSSM with large tan /3. Studies also show that 
the running behaviours can be quite different for Dirac or Majorana neutrinos i7|] , and the 
additional Majorana phases may have intrinsic behaviour in the evolution |20|]. In this 
section, we proceed to consider the one- loop RGEs of |V^jP, ^1' and for both the Dirac 
and the Majorana neutrinos. Note that all these parameters are rephasing invariant and 
independent of any particular parametrization of the MNS matrix V. Distinguishable RGE 
running effects between Dirac neutrinos and Majorana neutrinos are discussed in detail. 



A. Dirac Neutrinos 



If neutrinos are Dirac particles, their Yukawa coupling matrix must be extremely 
suppressed in magnitude to reproduce the light neutrino masses of 0{1) eV or smaller at 
low energy scales. The running of Y^^ from the electroweak energy scale A-^-^ to a superhigh 



energy scale A is governed by the one-loop RGE 



2l| 



167r' 



dt 



2aj,oo + C [{YiY,^) CO + CO (y,Y,^ 
10 



(29) 



where ui = Y^Xv is a Hermitian quantity, i = In (///A) with // being an arbitrary renormahza- 
tion scale between and A, Y"; is the charged-lepton Yukawa couphng matrix, C — —1.5 
(SM) or C = 1 (MSSM) and ~ -0.45^f _ 2.25^| + Sy^ (SM) or ~ -0.6c/f - 3^1 + Sy^ 
(MSSM). Here (^^ and g2 are the gauge couphngs, stands for the top-quark Yukawa cou- 
phng. In writing out Eq. (29), we have safely neglected those tiny terms of (9(a;^). 

Without loss of generality, we choose the flavor basis where is diagonal and real (pos- 
itive): Yi — diag l^/g, y^, y^}. In this basis ui can be diagonalized by the unitary transfor- 
mation V'^uV — Cj = diag{y^, y|, y|}, where V is just the MNS matrix and at A^^ the 
Dirac neutrino masses are m,i = vy^ (SM) or = vy^ sin /S (MSSM) with v ^ 174 GeV. One 
may use Eq. (29) to derive the explicit RGEs of neutrino masses and the MNS matrix V. 
By taking some lengthy but not complicated calculations we can further derive the RGEs of 
those rephasing invariant parameters with the help of Eqs.(2), (5) and (10). Here we simply 
give the resulting one-loop RGEs of \V^^\'^, $ and J. During the derivation, we have taken 
the approximation of r-lepton dominance. In other words, the contributions of y^ and y^ to 
all of the RGEs are neghgibly small and are safely neghgted. 
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(32) 



Since we have Arn^i ^ ^^31 ~ ^^32? "we can infer from Eq. (30) that in the standard 
parametrization 6^2 = arctan(| V^jl/lKil) is in general more sensitive to the radiative correc- 
tion than the other two mixing angles 9^^ = arcsin(| VggI) and 6^33 = arctan(| V^gl/I V^gl). This 
result is also true if the neutrinos are Majorana particles which can be easily seen from Eq. 
(41) and is consist with the analyses in many previous papers 



Since the absolute mass scale of three neutrinos and the sign of Amg^^ remain unknown, 
we further consider three typical patterns of the neutrino mass spectrum: normal hierarchy 
(NH), inverted hierarchy (IH) and near degeneracy (ND). 
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Normal Hierarchy mg ^ ^ 171^ ~ 0, ~ ^f^m^^ and mg y^A 
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In this neutrino masses limit the one-loop RGEs of $ and J can be approximately 
expressed as 
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(34) 



dt ^ -~-T ^ v,-rl, ,-r3l 

where in the next leading order terms lead by IS.vn^^l IS.vn^^^ we preserve only those 
terms inversely proportional to | V^g | . Taking into account of the smallness of | V^g | , the 
contribution of these terms may be comparable with the leading order terms. 

Some discussions are in order. 

— The angle is most insensitive to the radiative correction in the leading order 
if neutrino masses are of normal hierarchy. 

— In case of the MSSM with large tan/3, Amg^/Amg^ increases significantly with 
the increase of energy scale [2l|], therefore $^2) *^t1 ^'^'^ *^r2 probably 
having significant evolutions, especially when |V"gg| takes a small value. While in 
the SM, all nine angles receive only small radiative corrections. 

— The one-loop RGE of J is proportional to J itself, therefore its evolutions are 
apparently opposite for positive and negative J . 

— Consider the evolution of J from to a superhigh energy scale. We can find 
from Eq. (12) that \V-,x^ - jV^^gP is negative at In the SM, C < 0, we 
can expect that \ J\ only slightly increases during the evolution. In the MSSM, 
C > 0, and we can find from Eq. (30) that |V^P increases with the evolution 
while |V^p decreases. Therefore |V^iP — |V^3p keeps negative which indicates 
that \J\ will go approaching zero during the evolution. 
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Inverted Hierarchy -y/— Am|i m2 ~ » f=:i 

In this neutrino masses limit the one-loop RGEs of $ and J can be approximately 
expressed as 
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Note that, in the case of IH, Am§i is negative. 



• Near Degeneracy ~ m2 ~ m^^ and Am|2 ~ Am|;^. 

If three neutrino masses are nearly degenerate, the one-loop RGEs of $ and J can be 
approximately expressed as 
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(38) 



We can easily find that the RGEs of and J7 in the cases of IH and ND are ahke if 
neutrinos are Dirac particles. Eqs. (35) and (36) can be obtained from Eqs. (37) and (38) 
by simply choosing ml = —Aml-^. We can expect from above equations (35) - (38) that all 
nine angles and the Jarlskog J7 may have large evolutions, especially in the MSSM with 
large tan /3 if the neutrino mass spectrum is ND or IH. 



B. Majorana Neutrinos 



Majorana neutrino masses are believed to be attributed to some physics at a superhigh 
energy scale A, e.g., the seesaw mechanisms. But all these new physics point to the unique 
dimension-5 Weinberg operator in an effective theory after the corresponding heavy particles 



are integrated out 
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or 



C 
C 



1 ^ 



l^H -K- H^ll + h.c. , 



^ l^H, ■ K ■ Hjll + h.c. 



(SM) 
(MSSM) 



(39) 



15 



which lead to the effective Majorana neutrino mass matrix = nv"^ (SM) or M^^ = 
^ sin^ /3 (MSSM), with tan /3 denotes the ratio of the vacuum expectation values of two 
MSSM Higgs doublets. Here A is a cut off energy scale stands for the energy scale of new 
physics. The evolution of n from A down to the electroweak scale A^-^ is formally inde- 
pendent of any details of the relevent model from which k is derived. Below A the energy 
dependent of the effective neutrino coupling matrix k, is described by 

T" 



, d/t 
dt 



IGtt^— = a^K + C 



(40) 



at the one-loop level [23|, where ^ —'^gi + ^Ut + (SM) or ^ — l-25'i — Qgi + Qyf 
(MSSM) with A denotes the Higgs self-coupling in the SM. 

Similarly, one may use Eq. (40) to derive the explicit RGEs for neutrino masses and MNS 
matrix in the flavor basis where 1^ is diagonal and real. In this basis, we have k = VkV^ 
with k = diag {k^, /tg, Kg} where V is just the MNS matrix and at A^-^^ Majorana neutrino 
masses are = v'^n^ (SM) or = f^KjSin^/3 (MSSM). Then we can further calculate 
the RGEs of |V^jp, ($) and J'. We give only the concise results here, where again the 
excellent approximation of r-dominance are taken. 
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^ ^.12 (i^.2r - i^.ir) - ^2 (iK2r - ikiD] 



'^'^121 



[1,.23 (I V^l' - \VM - V23 - } . (43) 

Prom Eq. (20), we have 

where /3 and 7 are the next two flavor indices right after a. By using this equation, the 
one-loop RGE of ^-matrix for the Majorana neutrinos can then be easily obtained from Eq. 
(42). 

We can clearly see that the parameters '^a^ij and ^apij which are associated with the 
Majorana phases in V and not related to the one-loop RGBs for Dirac neutrinos are involved 
in the RGEs for Majorana neutrinos. These terms could dominate over others and determine 
the running behaviours of |V^j|, and J if the Majorana phases are properly chosen. It 
is well known that for Dirac neutrinos, if J is zero at some scale, it will keep vanished at 
any energy scale. However, we can see from Eqs. (42) and (43), for Majorana neutrinos, 
even if JT" = (no Dirac type CP violation) at some energy scale, $„j and J can still receive 
significant radiative corrections only if not all the ^^.^ are zero (i.e., Majorana type CP 
violation exists). 

Here we give the approximate RGEs of ^^.^ and J in three limits of neutrino mass 
hierarchy: NH, IH and ND. 
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• Normal Hierarchy 771,3 » 777.2 » ~ 0, m2 ~ V Attt.^^ and 7773 fa ^/^m^ 

In this neutrino masses limit the one-loop RGEs of ^ and J can be approximately 
expressed as 
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(45) 



+ 2Cy: 
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(46) 



Again, for terms lead by I\vr?2\l ^"^^xi reserve only those terms inversely propor- 
tional to |V^3|. By using Eq. (44), the one- loop RGEs of in the limit of NH can 
be easily derived from Eq. (45): 
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(47) 



Here the symbol x stands for terms that are not inversely proportional to |V'e3|, which 
lead only mild corrections to the angle matrix in the case of NH. 
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Comparing Eqs. (47) and (46) with Eqs. (33) and (34), we can find that the terms 
lead by \J Am|]^/Am|]^ may lead to very different running behaviours of Major ana 
neutrinos compared to the Dirac neutrinos. For some specific pattern of V ^ terms 
with '^ctfiij can be dominating and even change the evolution directions of and J . 



Inverted Hierarchy ^J^^Krr^^ 777.3 ~ » mg f=:i 

In this neutrino masses limit the one-loop RGEs of $ and J can be approximately 
expressed as 
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dt 
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+ ^(3iK3r-i)} ■ 

Note that, in the case of IH, A777I1 is negative. 



(49) 



Near Degeneracy rrig ~ ~ and Am|2 ~ ^^31 ■ 

If three neutrino masses are nearly degenerate the one-loop RGEs of and J can be 
approximately expressed as 
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Eqs. (45) - (51) indicate that the RGE running behaviours of ^ and J for Majorana 
neutrinos are very different to that for Dirac neutrinos (see Eqs. (33) - (38)). We can 
find that three neutrino masses are nearly degenerate, the running behaviours of $, \I' and 
J depend on the interplay of several terms (see Eqs. (37), (38), (50) and (51)) and are 
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very sensitive to the sign of Am|i no matter whether the neutrinos are Dirac or Majorana 
particles. 

IV. NUMERICAL ANALYSIS AND DISCUSSION 

The running behaviours of the above mentioned rephasing-invariant quantities arc numer- 
ically illustrated by assuming A ~ 10^"^ GcV, which is the typical scale of the conventional 
seesaw mechanisms and is very close to the scale of the grand unified theories. We chose 
several sets of typical values of the angle matrix elements in Eq. (9) (Eq. (16)) which are 
allowed by current 3a experimental data at and calculated the RGE running effects of 
the Dirac (Majorana) angle matrix $ (^) and the Jarlskog invariant JT" from A^-^^ up to A. 
For each set of inputs, we consider four typical pattterns of neutrino mass spectrum: i) NH 
(m^ ~ 0) , n) IH (rrig ~ 0) , iii) ND with Amli > and iv) ND with Am|]^ < 0. In our 
numerical calculation, Am|i = 7.59 x 10~^ eV^ and Amg^ = ±2.4 x 10"^ eV^ have been 
taken as the typical inputs at A^^ and in case iii) and iv) we have chosen m^^ = 0.2 eV. 
We carry out our numerical calculation in the framework of either the SM or the MSSM, 
where the Higgs mass = 140 GeV in the SM and the parameter tan ;5 = 10 or 50 in the 
MSSM have typically been input. Our numerical results and the corresponding inputs are 
summarized in Tables I - IV and Figs. 2-5. 

In Table I-III, for the same set of inputs, we calculated the radiative corrections to J", $ 
and at A in the SM (Table I), the MSSM with tan ^ = 10 (Table II) and the MSSM with 
tan ^ = 50 (Table III) respectively. If the neutrino mass spectrum is NH or IH, angles of 
$ and ^ can receive non-negligible radiative corrections (larger than 1°) only in the case of 
MSSM with tan/3 = 50. If three neutrino masses are nearly degenerate, all the angles of $ 
and \I' can receive significant radiative corrections especially in the MSSM and we can find 
that their running behaviours are very sensitive to the sign of Amli for either Majorana or 
Dirac neutrinos. 

The running behaviours of the Jarlskog invariant J in the SM and the MSSM with 
tan /3 — 10 are shown in Fig. 2 and 3 respectively for both the Dirac and Majorana neutrinos. 

The Jarlskog J' can receive significant radiative correction if three neutrino masses are nearly 
degenerate. As already mentioned in Section 111, if neutrinos are Majorana particles, the 
running of is very sensitive to the sign of Am|]^ in case of ND. For this specific set of 
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inputs (with negative ^7), we can see that if A.m^i < 0, will decrease when running from 
Ag-y^ to A. If A.ml^ > 0, J will increase and in case of MSSM with tan /3 — 50, J can even 
run above zero and evolve to a positive value (see Case (a) in Fig. 5). 

In the framework of MSSM with tan /3 = 50, we consider another two sets of inputs 
with J' > and = respectively. The corresponding results are shown in Table IV 
and V. Fig. 4 and 5 illustrate the evolution of J" in these three special cases: a) J" < 
(Table III), h) J > (Table IV) and c) J = (Table V) for Dirac and Majorana neutrinos 
respectively. Note that the input values of J and ^^.^ at Ag^ in Table IV have the same 
absolute values but opposite signs as that in Table III. If neutrinos are Dirac particles, the 
RGE running behaviours of J' and in these two cases are entirely opposite. We can 
find that at A ~ 10^'' GeV, J' and <I>^j in Table III and IV still have the same absolute 
values but opposite signs. For Dirac neutrinos, \J'\ is always decreasing and approaching 
zero in the MSSM when running from A^^ to A. It means that when evolving from A^^ 
to A, positive J decreases while negative J increases on the contrary, but the signs of J 
are not changed. However, if neutrinos are Majorana particles, the running behaviours of 
the Dirac- type CP- violating parameters JT" and strongly depend on the Majorana- type 
CP- violating parameters A point need to be pointed out is that the running behaviours 
of i7 in a), b) and c) three cases are somehow similar as shown in Fig. 5. This is because 
of these three sets of inputs correspond to the same values of the Majorana angles p and a 
in the Standard Parametrization. By choosing a different set of inputs of "if^^, we are able 
to make a very different evolution of J for the same set of ^^j. A more special case is that 
J = but the Majorana CP violation is nonzero as shown in Table V. If JT" = is input at 
^EW' ^ ^^^^ keep vanishing in case of Dirac neutrinos. But for Majorana neutrinos, can 
evolve to a nonzero value at A, which indicates that the UTs expand from a line and the 
Dirac-type CP violation is radiatively generated from the Majorana-type CP violations. 

As we have declaimed in the beginning section, both the Dirac-type and the Majorana- 
type CP violation can be illustrated by the six UTs. Therefore, the RGE running of the 
Dirac angle matrix $ and the Majorana angle matrix ^ also correspond to the evolutions 
and the rotations of the UTs in the complex plane. We choose two examples in the MSSM 
with tan fS = 50 where the Majorana neutrino masses are nearly degenerate and show in Fig. 
6 {Amli > 0) and 7 (Amg^ < 0) how the six UTs evolve from A^^ to A, where triangles 
with thicker sides are at higher energy scale. The common area of the six UTs equals to 
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i7/2. Since the orientations of A^, and A^ have no physical meaning, we simply choose 
one side of each triangle {V^^^*^, K3K3 V^^V*^) to lie on the x-axial and point to the 
origin. In this way, the triangle lie above the x-axial is clockwise and corresponds to a 
positive and the triangle lie below x-axial is anti-clockwise corresponds to a negative 
on the contrary. Then we can clearly see from Fig. 6, when evolve from to the high 
energy scale A, the Jarlskog changed its sign as illustrated in Table III. 

V. SUMMARY 

In summary, we introduced the concepts of the Dirac angle matrix $ and the Majorana 
angle matrix ^ for Dirac and Majorana neutrinos respectively and show that the angle 
matrix carries equivalent information to the complex mixing matrix itself, but with the added 
advantage of being basis and phase convention independent. Our prescription works for any 
number of fermion generation. We further calculated the one-loop RGEs of ^ and some 
other rephasing invariant parameters. Numerical analyses are carried out for illustration. 
We find that apparently different from the case of Dirac neutrinos, for Majorana neutrinos 
the RG-evolutions of ^ and J strongly depend on the Majorana-type CP-violating 
parameters and are quite sensitive to the sign of Am|]^. They may receive significant radiative 
corrections in the MSSM with large tan f3 if three neutrino masses are nearly degenerate. 

Of course, the numerical examples presented in this work are mainly for the purpose of 
illustration. The point is that the nature and the mass spectrum of neutrinos determine 
the RGE running behaviours of those rephasing invariant parameters which may be crucial 
for building a realistic neutrino model. Our analysis complement those previous studies of 
radiative corrections to the physical parameters of Dirac and Majorana neutrinos and are 
helpful for building a realistic neutrino mass model at a high energy scale. 
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Appendix A: Reconstruct the Leptonic Mixing Matrix V from the Dirac Angle 
Matrix $ or the Majorana Angle Matrix 

In Section II we have shown that the Dirac (Majorana) angle matrix carries equivalent 
independent real parameters as that of the MNS matrix V for Dirac (Majorana) neutrinos. 
We now further prove the equivalence of the mixing matrix and the angle matrix in the 
presence of the Dirac- type CP violation (which indicates that none of the moduli \V^^\ is 
zero) by presenting the way of re-obtaining the mixing matrix starting from the angle matrix 
in case of three generation neutrinos. We shall take it as given that the matrix of mixing 
moduli is essentially equivalent to the complex mixing matrix {2^, and content ourselves in 
the first instance with showing how to obtain the mixing-matrix moduli \V^^\, starting from 
the angles. 



1. Dirac Angle Matrix $ 



The 3x3 Dirac angle matrix $ is defined by Eqs. (9) and (10), and its nine matrix 
elements satisfy the normalization conditions of Eq. (11). The Dirac angle matrix can be 
manipulated to yield the magnitudes of mixing matrix elements. Firstly, we introduce the 
sin-matrix 
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and then define certain products of sines Sq,,: 



sin sin sin sin , 



(A2) 



multiplying together the four sin $ entries in the same row and column as sin excluding 
sin$j^^ itself. Clearly every mixing modulus-squared except |V^jP enters in the denominator 
of the product S^.^, whereby the S^.^ must be proportional to |V^jP, 

1 



\V • 



(A3) 
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The relevent normalising factor N may be obtained by summing over any row or column 
(or indeed over both rows and columns) 



Then we have 



As for the Jarlskog J ^ we can find from Eq. (Al) that 
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together with Eq. (A4) and (A5), we can obtain 
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J = ^^^^ ■ (A7) 



Then we conclude that the ^-matrix is equivalent to the complex mixing matrix V if neu- 
trinos are Dirac particles. 



2. Majorana Angle Matrix ^ 

The Majorana angle matrix is defined as 
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from which we can easily obtain the $-matrix defined in Eqs. (9) and (10) by using the 
relation of Eq. (20). Follow the same procedure as that in above section, we can obtain the 
|y|-matrix and the Jarlskog J from the matrix. Two extra Majorana phases can also be 
deduced from the ^-matrix itself. 
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For a specific parametrization of V, for example the Standard Parametrization 
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three mixing angles and the Dirac phase 8 can be determined from the -matrix and the 
Jarlskog J: 
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Both the absolute value and the quadrant of 5 can be determined by Eqs. (A13) and (A14). 
Two Majorana phases p and a can also be easily obtained through 



P = <^ - ^e2 , 



(A15) 
(A16) 
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FIG. 1: Schematic diagrams for six leptonic unitarity triangles in the complex plane, where each 
triangle is named by the index that does not manifest in its three sides. 
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TABLE I: Radiative corrections to the Jarlskog J ^ the Dirac angle matrix $ and the Majorana 
angle matrix at A ~ 10^^ GeV in the SM (with = 140 GeV) where we choose — 34°, 
6I23 = 46°, 6*13 = 7°, (5 = -90°, p = 120° and a = 60° (in the standard parametrization of V) as 
typical inputs at the electroweak energy scale A^^. 
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-9.40° 


-9.42° 
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— 2U.5U 


-20.50° 


-20.49° 


-20.48° 


-20.30° 
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-20.47° 


-20.15° 


-20.10° 


/Fx 
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-150.10° 


-150.10° 


-150.02° 


-150.29° 


-150.10° 


-150.11° 


-150.27° 


-150.32° 




-00.40 


-85.46° 


-85.43° 


-84.93° 


-84.91° 


-85.46° 


-85.47° 


-85.64° 


-85.64° 




on 1 no 
— oU.lO 
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TABLE II: Radiative corrections to the Jarlskog J ^ the Dirac angle matrix $ and the Majorana 
angle matrix at A ~ 10^"^ GeV in the MSSM (with tan/3 = 10) where we choose Qy^ = 34°, 
6»23 = 46°, 6li3 = 7°, (5 = -90°, p = 120° and a = 60° (in the standard parametrization of V) as 
typical inputs at the electroweak energy scale A^^. 
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TABLE III: Radiative corrections to the Jarlskog J ^ the Dirac angle matrix $ and the Majorana 
angle matrix at A ~ 10^"^ GeV in the MSSM (with tan/3 = 50) where we choose Qy^ = 34°, 
6»23 = 46°, 6li3 = 7°, (5 = -90°, p = 120° and a = 60° (in the standard parametrization of V) as 
typical inputs at the electroweak energy scale A^^. 
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-118.75° 
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-00.40 
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-109.98° 


16.84° 


-60.88° 


-85.24° 
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-14.67° 
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-81.18° 


-56.79° 


129.41° 


-100.27° 
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-83.31° 


-77.09° 


-78.01° 




— 14.44 


-14.44° 


-13.23° 


33.75° 


-18.85° 


-15.09° 


-35.03° 
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-87.32° 
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— oo.i4 
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-11.95° 


70 QdO 
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-28.75° 


-45.50° 


100.87° 
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-29.09° 


-25.58° 
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-50.30° 










*e3 
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57.84° 


71.08° 


101.04° 


78.63° 










Ml 


64.86° 


64.15° 


71.85° 


79.48° 


128.46° 










*;.2 


70.60° 


71.77° 


55.09° 


-6.43° 


117.75° 










*m3 


-135.46° 


-135.92° 


-126.94° 


-73.05° 


113.79° 










*rl 


-124.54° 


-124.37° 


-115.52° 


-95.97° 


-147.45° 










*r2 


-129.90° 


-127.91° 


-118.79° 


-151.32° 


-130.03° 










*r3 


-105.56° 


-107.72° 


-95.69° 


-112.71° 


-82.52° 
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TABLE IV: Radiative corrections to the Jarlskog J ^ the Dirac angle matrix $ and the Majorana 
angle matrix * at A ~ 10^"^ GcV in the MSSM (with tan/3 = 50) where we choose Qy^ = 34°, 
6»23 = 46°, = 7°, 6 = 90°, p = 120° and a = 60° (in the standard parametrization of V) as 
typical inputs at the electroweak energy scale A^.^. 







A ~ lO^'' GeV 


Majorana Neutrinos 


Dirac Neutrinos 


NH 

m^^ ~ 


IH 

mg ~ 


Near Degeneracy 


NH 

m^^ ~ 


IH 

mg ~ 


Near Degeneracy 


Amgi > 


Am^i < 


A 9 ^ 

Am§i > 


Am^i < 


J 


0.027812 


0.027568 


0.015023 


0.014327 


-0.027134 


0.026677 


0.010653 


0.000337 


0.000332 


%1 


9.40° 


9.13° 


4.46° 


4.94° 


-40.79° 


8.84° 


2.61° 


0.43° 


0.43° 




20.50° 


21.22° 


17.53° 


47.79° 


-125.26° 


20.48° 


58.64° 


88.07° 


101.29° 




150.10° 


149.65° 


158.01° 


127.27° 


-13.95° 


150.68° 


118.75° 


91.50° 


78.28° 


^.1 


85.46° 


83.82° 


35.79° 


11.03° 


-124.48° 


85.24° 


61.66° 


15.17° 


14.67° 




80.10° 


80.66° 


132.66° 


122.51° 


-43.25° 


79.67° 


83.31° 


77.09° 


78.01° 




14.44° 


15.52° 


11.55° 


46.46° 


-12.27° 


15.09° 


35.03° 


87.74° 


87.32° 




85.14° 


87.05° 


139.75° 


164.03° 


-14.73° 


85.92° 


115.73° 


164.40° 


164.90° 




79.40° 


78.12° 


29.81° 


9.70° 


-11.49° 


79.85° 


38.05° 


14.84° 


0.70° 




15.46° 


14.83° 


10.44° 


6.27° 


-153.78° 


14.23° 


26.22° 


0.76° 


14.40° 




150° 


148.49° 


111.98° 


94.10° 


-66.18° 












150° 


151.76° 


161.04° 


152.04° 


-20.50° 












60° 


59.75° 


86.98° 


113.86° 


86.68° 










* 1 


55.14° 


55.54° 


71.73° 


78.13° 


99.09° 












49.40° 


49.88° 


10.85° 


-18.26° 


148.01° 












-104.54° 


-105.42° 


-82.58° 


-59.87° 


112.90° 










*rl 


-115.46° 


-115.33° 


-103.81° 


-96.93° 


-121.70° 












-110.10° 


-108.90° 


-151.02° 


-150.-17° 


-157.25° 










*r3 


-134.44° 


-135.77° 


-104.57° 


-112.60° 


-81.05° 
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TABLE V: Radiative corrections to the Jarlskog J ^ the Dirac angle matrix $ and the Majorana 
angle matrix at A ~ 10^"^ GeV in the MSSM (with tan/3 = 50) where we choose Qy^ = 34°, 
6»23 = 46°, 6*13 = 7°, (5 = 0°, /9 = 120° and a = 60° (in the standard parametrization of V) as 
typical inputs at the electroweak energy scale A^^. 













A ~ 10^4 


GeV 












Majorana Neutrinos 


Dirac Neutrinos 






NH 


IH 


Near De 


generacy 


NH 


IH 


Near De 


generacy 






m^^ ~ 


777,3 ~ 


Am|i > 


Am|i < 


m^^ ~ 


777,3 ~ 


Amli > 


Amli < 


J 





0.001043 


-0.011926 


0.013144 


-0.036288 
















0° 


0.35° 


-3.67° 


4.56° 


-67.38° 


0° 


0° 


0° 


0° 




0° 


0.83° 


-12.41° 


40.09° 


-94.47° 


0° 


0° 


0° 


0° 




180° 


178.82° 


-163.92° 


135.35° 


-18.15° 


180° 


180° 


180° 


180° 




0° 


2.03° 


-26.13° 


12.83° 


-88.40° 


0° 


0° 


0° 


0° 




180° 


177.31° 


-145.45° 


128.01° 


-72.50° 


180° 


180° 


180° 


180° 




0° 


0.66° 


-8.42° 


39.16° 


-19.10° 


0° 


0° 


0° 


0° 




180° 


177.62° 


-150.20° 


162.61° 


-24.22° 


180° 


180° 


180° 


180° 




0° 


1.86° 


-22.14° 


11.90° 


-13.03° 


0° 


0° 


0° 


0° 




0° 


0.52° 


-7.66° 


5.49° 


-142.75° 


0° 


0° 


0° 


0° 




60° 


62.12° 


42.14° 


93.61° 


-35.78° 












-120° 


-121.42° 


-115.33° 


-12.72° 


-44.94° 












60° 


59.30° 


73.19° 


-80.89° 


80.72° 










* 1 


60° 


59.74° 


71.94° 


76.22° 


120.00° 












60° 


60.44° 


42.53° 


179.18° 


122.03° 












-120° 


-120.18° 


-114.47° 


104.60° 


117.97° 










*rl 


-120° 


-119.91° 


-111.73° 


-99.22° 


-127.38° 












-120° 


-118.73° 


-1-19.88° 


39.27° 


-152.11° 










*r3 


-120° 


-121.36° 


-98.39° 


59.95° 


-80.18° 
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FIG. 2: Running behaviours of the Jarlskog invariant J from to A ~ 10 GeV in the SM for 
both the Dirac and the Majorana neutrinos, where the corresponding inputs at A^^ and outputs 
at A can be found in TABLE I. In this figure, the sohd hncs stand for the case of NH with ~ 0, 
the dotted hues for the case of IH with 777,3 ~ 0' dashed lines for the ND case with Amg^^ > 
and the dash dotted hues for the ND case with Amg^ < 0. Two sohd hues in this figure are almost 
coincide. 
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FIG. 3: Running behaviours of the Jarlskog invariant J from to A ~ 10^^ GeV in the MSSM 
(with tan /3 = 10) for both the Dirac and the Majorana neutrinos, where the corresponding inputs 
at Ag^ and outputs at A can be found in TABLE II. In this figure, the sohd hues stand for the 
case of NH with « 0, the dotted lines for the case of IH with mg k, 0, the dashed lines for the 
ND case with Amg^^ > and the dash dotted lines for the ND case with Amg^ < 0. Two solid 
lines in this figure are almost coincide. 
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FIG. 4: Running behaviours of the Jarlskog invariant J from A^-^^ to A ~ 10^^ GeV in the MSSM 
(with tan /? = 50) for Dirac neutrinos, where the corresponding inputs at A^^ and outputs at A 
for case a), b) and c) can be found in TABLE III, IV and V respectively. The solid lines stand for 
the case of NH with m^^ ~ 0, the dotted lines for the case of IH with rrig ~ 0, the dashed lines for 
the ND case with Amg^^ > and the dash dotted lines for the ND case with Amg^ < 0. In this 
figure, the dashed lines and the dash doted lines are almost coincide in case I and II, and in case 
III all four lines are almost coincide. 
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FIG. 5: Running behaviours of the Jarlskog invariant J from A^-^^ to A ~ 10 GeV in the MSSM 
(with tan /3 = 50) for Majorana neutrinos, where the corresponding inputs at A^-^ and outputs 
at A for case a), b) and c) can be found in TABLE III, IV and V respectively. In this figure, the 
solid hues stand for the case of NH with ~ 0, the dotted lines for the case of IH with mg ~ 0, 
the dashed lines for the ND case with Arn-g^^ > and the dash dotted lines for the ND case with 
Amoi < 0. 
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FIG. 6: RG-evolutions of all six leptonic unitarity triangles in the complex plane from to 
A ~ 10^^ GeV in the MSSM (with tan^ = 50) for Majorana neutrinos, where triangles with 

thicker sides are at higher energy scale. Here the neutrino mass spectrum is ND with Ani^i > 
and the corresponding inputs at Ag^ and outputs at A can be found in TABLE III. 
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FIG. 7: RG-evolutions of all six leptonic unitarity triangles in the complex plane from to 
A ~ 10^^ GeV in the MSSM (with tan^ = 50) for Majorana neutrinos, where triangles with 

thicker sides are at higher energy scale. Here the neutrino mass spectrum is ND with Ani^i < 
and the corresponding inputs at Ag^ and outputs at A can be found in TABLE III. 
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